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Abstract 

We study theoretically the non-linear response properties of glass formers. We establish several 
general results which, together with the assumption of Time- Temperature Superposition, lead to 
a relation between the non-linear response and the derivative of the linear response with respect 
to temperature. Using results from Mode-Coupling Theory (MCT) and scaling arguments valid 
close to the glass transition, we obtain the frequency and temperature dependence of the non- 
linear response in the a and /3-regimes. Our results demonstrate that supercooled liquids are 
characterized by responses to external perturbations that become increasingly non-linear as the 
glass transition is approached. These results are extended to the case of inhomogeneous perturbing 
fields. 
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I. INTRODUCTION AND MOTIVATIONS 



Physical systems become very susceptible to external perturbations close to a phase 
transition. In several theoretical approaches the huge increase of the relaxation time of 
super-cooled liquid is traced back to the proximity to a phase transition 4l|. Therefore, it 
would be natural to expect that responses to external fields become singular also approaching 
the glass transition. However, standard linear responses (as well as correlation functions) 
are known to remain quite featureless. In this work, following up 1|, we show using general 
arguments that instead non-linear responses do increase. Furthermore, we work out detailed 
predictions within the Mode Coupling Theory of glasses and we show what is the relation 
between the growth of non-linear responses and the one of previously introduced probes of 
dynamic correlations. 

In order to grasp why only non-linear responses may grow approaching the glass transi- 
tion, it is instructive to recall the situation of standard critical phenomena and contrast it 
to the one of spin glasses. In standard critical phenomena there is a spontaneous symme- 
try breaking toward an ordered phase. The main consequence of the broken symmetry is 
that there are different but thermodynamically equivalent states in which the system can 
be found: they are all related by the symmetry at hand. For example, below the ferromag- 
netic transition the state with a given magnetization M and the one characterized by the 
magnetization — M are equivalent and the system will be in one or the other depending on 
boundary conditions, residual fields, etc. Switching on an external field that couples to the 
order parameter allows one to select one given state, e.g. a positive magnetic field selects 
the state with positive magnetization in ferromagnetic systems. As a consequence, it is easy 
to understand that just before the transition, where the system is at the brink of developing 
long range order, the response to an external field is huge and actually becomes infinite at 
the critical point. The situation is more subtle for systems characterized by the appearance 
and growth of amorphous long range order, as it is the case in spin glasses and is conjectured 
to be the case for glasses 

aa a. 

The only diverging linear susceptibility would be the 
staggered one: the response obtained by using a small field with spatial modulation corre- 
lated to the one of the amorphous state. Of course, this is not doable in experiments and 
not even in simulations since it is impossible to guess the external field which would impose 
a given amorphous long-range ordered state because it is amorphous as the state it pins (an 
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external field with an uncorrelated spatial modulation, e.g. an homogenous field, does not 
allow to select one state from the others and so it does not lead to singular responses). For 
spin glasses, the square of the staggered linear response can nevertheless be probed, since 
it is related to the static third order non-linear susceptibility 0, [(], [7] , which was shown 
to diverge at the transition. Similar arguments [l| (see also below) adapted to supercooled 
liquids - where the disorder is not quenched but self-induced - suggest that linear responses 
will be featureless but non-linear susceptibilities should instead increase approaching the 
glass transition. The main difference between spin glasses and structural glasses is that in 
the latter case one has to focus on dynamical non-linear responses. It was argued in [l| , 
on the basis of physical and heuristic arguments, that the non-linear dielectric susceptibility 
X3(a>) should exhibit a growing peak around ur a = 1, while Xzip = 0) should remain trivial, 
in contrast with the case of spin-glasses (r a is the standard notation for the relaxation time 
of super-cooled liquids). However, the detailed shape of X3( u ) i n the glassy region is beyond 
the grasp of those heuristic arguments. Since the corresponding experiments are currently 
being performed a, ISI], it is quite important to get more precise predictions on the expected 
shape of ^3(0;). This is the primary aim of the present study, where we obtain for the 
first time, using general arguments made more precise in the context of the Mode Coupling 
Theory (MCT) lOj of the glass transition, some precise information on the non-linear sus- 
ceptibility, in particular concerning both its frequency and temperature dependence. Our 
results can be generalized to other non-linear responses such as non-linear compressibility 
and non-linear rheological responses. These should be measurable for colloids approaching 
the glass transition. 

It is important to relate the results on non-linear responses to the ones obtained re- 
cently on dynamical correlations. This is the second aim of our work. For purpose of self- 
consistency, we recall below some definitions and results that we will use in the following 
sections. 

It was established in the last decade that glass- forming liquids become more and more 



dynamically correlated approaching the glass transition ll|. This phenomenon, related to 



an increasing heterogeneity in the dynamics, is remarkable since, for the first time, some 
type of spatial correlation has been clearly connected to the slowing down of the dynamics 
of supercooled liquids. In order to unveil the existence of dynamical correlations, one has 
to focus on some local probe of the relaxation dynamics - typically a two point function 
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0(x,t) = (f>(x,t)(f)(x,Q) where 
and time t, or the mobility field 



x,t) could be the density fluctuation, 6p(x,t), at position x 



121 ] . A measure of dynamical spatial correlations is obtained 

nn 



by considering the four point correlator G 4 (x — y;t) = (0(x,t)0(y,t)) c [13|, |14J, where the 
brakets denotes the usual connected average over dynamical histories. By analogy with 
standard critical phenomena, a 'susceptibility' 42J may be naturally defined by spatially in- 
tegrating G4: Xi(t) = I dxG&{x,t). This 'dynamical susceptibilit y', ^ 4, has been intensi vely 
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3,Q, 
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2l|. 



studied in the past few years, both theoretically and numerically 
One finds that Xa reaches a peak value for time scales of the order of the relaxation time 
of the system r a , and the height of this peak increases as the temperature is reduced, as 
a clear sign of the growth of some dynamical correlation length as the glass transition is 
approached. From an experimental point of view, however, four-point correlation functions 
are very difficult to measure directly, except in cases where one can monitor the trajectory of 
individual particles - for example granular and colloidal systems where can be measured 
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directly and again shows interesting features as the system jams [22, 

Another interesting quantity, extensively studied in the context of dynamic heterogene- 
ity in the past few years, is the derivative of the standard two-body correlation C[f) (or 
susceptibility xi{w)) with respect to the temperature (or the density) - a quantity called 
Xt = TdC(r)/dT (or xp = p9C(r)/dp) in 26|. This is a non-standard linear dynamical 
response and is clearly an easily accessible quantity, which also shows a peak at times of the 
order of r a and whose height grows as the temperature is lowered 26|, |27[ . This has lead to 
direct estimates of the size of the dynamical correlation length in supercooled liquids and 



glasses 26|, |27|, |28( . The relation between Xi an d Xt is however highly non-trivial and has 



been investigated thoroughly in [19|, [20(. It was realized in these papers that the existence of 
conserved quantities (energy, density) crucially affects the properties of xa-i which depends 
both on the thermodynamic ensemble (NVE vs. NPT for example) and on the dynamics 
(Brownian vs. Newtonian for example). The true glassy correlation length, on the other 
hand, does not depend on these choices, and therefore the direct interpretation of Xi(t) i n 
terms of a correlation volume is somewhat obscured. At a deeper level, the basic ingredient 
leading to the critical behaviour of Xi turns out to be entirely contained in the response 
function xt itself, as the field theoretical analysis of 19| explicitely demonstrates and the 
numerical results presented there fully confirm. For example, for Brownian dynamics or for 
Newtonian dynamics in the NVE ensemble, Xi ~ Xt, whereas for Newtonian dynamics in 
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the NVT ensemble, one rather finds %4 



sec 
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20| for a detailed discussion. 



A. Summary of results 

One of the two main results of this work consists in obtaining the precise relation between 
non-linear responses and probes of dynamical correlations introduced previously in the lit- 
erature. In particular, we shall focus on the third-order non-linear dielectric response 
and work out its relationship with the dielectric dynamical susceptibility xt- The latter is 
known to be related to X4> see 19j for a detailed discussion. We will establish simple identi- 
ties between \z an d Xt which hold whenever Time- Temperature Superposition (TTS) holds, 
i.e. whenever all dependence of the linear response of the system on external parameters 



(temperature, density, electric field, . . . ) comes through the dependence of the relaxation 

lin MCT, where a true dynamical 
10] . More precisely, we will show 



time on these parameters. This is a strict statement wit 
phase transition takes place at a finite temperature T c 
that at low enough frequency (in fact much smaller than the inverse of the relaxation time), 
the following relation holds: 

3(Xs(o;))«^Xr(2a;), (1) 

where k is a slow varying function of the temperature, see the following sections for a precise 
definition. This relation is actually a consequence of a more general one relating Xz{^) to 
the standard linear response Xi^X which is also valid only at low frequency and reads: 

XaH « k dT , (2) 

These simple relations do not extend over the whole frequency domain, however one could 
argue that they should do so in a scaling sense (as long as the TTS is valid). A general 
proof based on field theory along the lines of [jj]] is left for the future. In this work, we have 
checked that within MCT this is indeed what happens. Furthermore, within MCT, we have 
also worked out the complete critical behavior of x 3 , which is sketched in Fig. [TJ 

• In the a-regime, i.e. u ~ l/r Q ~ e 1 ^ 2a+1 ^ 2b /t , X3( u ) growths and reaches its max- 
imum, of height of order 1/e, after which it decreases as uJ~ b Tp at large uj. In this 
regime, one has the scaling form: Xzi^) — -G{uT a ). 

• At the crossover between the early a-regime and late /3-regime Xzd^) is of order l/\/e. 




log CO 



FIG. 1: Sketch of log |x3(w)| as a function of logo;, showing five different frequency regimes: 
ur a -C 1, uJT a ~ 1, Tp/r a <C uJTp <1 (e = T — T c ), urp 3> 1, ujtq ~ 1. Note that the low frequency 
limit is non zero but much smaller than the peak value for T close to T c . 

• In the /3-regime, i.e. tu ~ 1/rp ~ e 1 / 2 "/^, Xsl 1 ^) dicreases as oJ~ b Tp at small and as 
uj~ a Ti3 at large tu, with scaling form Xzi *) = ^i^i UJT l3) 

Exponents a and b are well known critical exponents of MCT, which characterize respectively 
how the correlators decay and exit of the plateau; tq is a microscopic relaxation time. All 
along this paper, e = T — T c , the distance from the Mode Coupling critical temperature T c , 
will be our control parameter. We remark that the existence of the peak and the decrease 
at low frequency is a non trivial prediction since it is in contrast to what happens for the 



(trivial) non-linear response of uncorrelated Brownian dipoles 29| and for spin glasses (the 



decrease with an exponent three at high frequency sketched in Fig. [T] is instead trivial and 
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will be discussed later). 

For symmetry reasons, the quadratic non- linear susceptibility Xii^) is zero f° r unpo- 
larized systems. This would not be the case for example for polarized systems, or when 
considering the response to a density perturbation. Indeed, the non-linear response to a 
density perturbation contains a quadratic term. The arguments presented below make it 
clear that in that case X2{u) itself is directly related to xt- We therefore expect that all the 
results presented below can be straightforwardly generalized to describe the long- wavelength 
(q — > 0) non-linear compressibility (or more general non-linear rheological responses) of su- 
percooled liquids and colloids. Moreover, by adaptin g th ese arguments to the explicit results 
obtained for Inhomogenous Mode Coupling Theory [30|] we will also extend these results to 
finite wavevectors q and find that whenever TTS holds, Xz{ u ) takes in the a-region the 
scaling form conjectured in 

X3 (uj)^C-"G^r a ), (3) 

where £ Q is the dynamical correlation length which also appears in Xt, and rj a certain 
critical exponent that within MCT is equal to minus 2. 



B. Organisation of the article 

The organisation of the paper is as follows. We first introduce the theoretical framework 
needed to deal with non-linear responses to an external field and establish some general 
relations between different quantities that naturally appear (Section [III)- We then exploit - 
when they exist - the Time- Temperature superposition (TTS) properties of the correlation 
function of glassy systems to establish in Section [TTT1 general relations, valid at low frequency, 
between X3 an d the temperature derivative of the linear response. The critical behavior of 
Xsi^) is obtained in Section [IV] using scaling arguments and within MCT. We finally dis- 
cuss the extension of these results to spatially inhomogeneous perturbing fields and beyond 
MCT. We end by a conclusion with open problems, possible extensions and experimental 
suggestions. 
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II. NON-LINEAR SUSCEPTIBILITY: GENERAL FRAMEWORK 



In this Section, we introduce the formalism needed to deal with non-linear response and 
establish a general relation between the non-linear susceptibility and a dynamical response 
function recently introduced in the literature, which was argued to capture the critical spatio- 
temporal correlations of the dynamics in the glassy region. In order to remain close to recent 
and ongoing experiments on glycerol, we use below the language of dielectric susceptibility. 
However, as mentioned in the Introduction, our arguments and results can be extended to 
more general non- linear susceptibilities (mechanical, magnetic, etc.). 



A. Linear and non linear response: small field expansion 

Let us consider a dipolar molecular liquid in presence of a small external electric field 
oscillating at frequency u in the ^-direction. We denote it as: 

E(t) = zE(t) = zE cos(wt), (4) 

where z is the unit vector in the z direction and E(t) = E cos(ut) is its ^-component with 
peak field amplitude E 

When the external field is sufficiently small, the polarization vector (per particle) P(t, E) 
can be expanded in powers of E. In the following we will denote P(t, E) its ^-component. 
Due to the rotational symmetry in the x-y plane the other components are identically zero. 
Furthermore, because of the up-down symmetry in the z direction, the polarization must be 
an odd function of E, i.e., P(t, —E) = —P(t,E). As a consequence, the expansion of P in 
powers of E contains only odd terms: 

P(t, E) = P x {t)E + P 3 (t)E 3 + 0(£ 5 ), (5) 

where Pi(t) and Ps(t) can be expressed as functional derivatives of the polarization with 
respect to the external field: 

8P(t) 



Pl(t) EE / d*i 

Ju<t 



tl <t SEih 



cos(a>£i) (6) 

E=0 



Ps(t) = \ I dfidfedV 



cos(a>ti) cos(a;t2) cos(a>£3) 

E=0 



It is important to remark that the linear and non-linear response kernels in the above 
integrals are time translation invariant (TTI), i.e. they do not change if all time variables 
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are shifted by the same amount. This comes from the fact that they are equilibrium response 
functions, measured in absence of the external field. Using this result and the specific form 
of the external field, Eq. (jlj), one finds: 

P(t, E)=E^ (xiM e™) + ^ (xxM e™ + * 3 (") e 3 ^) + 0(E 5 ). (7) 

which defines the usual frequency dependent linear susceptibility, Xi(^), an d the frequency 
dependent non-linear susceptibility, Xzij^), while Xi,2(^) is the E 2 correction to the first 
harmonic susceptibility Xiil^)- The non-linear susceptibility Xzi^), which is the quantity 
will focus on throughout this paper, can be accessed experimentaly by filtering 0(E 3 ) terms 
at frequency 3a;. 

Following the same procedure, one can expand in powers of the electric field the (z- 
component) polarization correlation and linear response functions of the system driven by 
the electric field E(t). The up-down symmetry in the z direction implies that they both are 
even functions of E. Therefore their expansion in power of E contains only even terms: 

C(t,t') = C (t,t')+C 2 (t,t')E 2 + 0(E 4 ) 

R(t,t') = R (t,t') + R 2 (t,t')E 2 + 0(E*). (8) 

Cq and i?o are the unperturbed correlation and response functions in absence of the external 
field. At equilibrium, they are functions only of the time difference r = t — t' > 0: Co{t, t') = 
Co(t — t') and Ro(t,t') = Ro(t — t'). Moreover, the Fluctuation-Dissipation theorem (FDT) 
holds for the unperturbed correlation and response functions: 

The second-order correlation and response functions appearing in Eq. (jSJ) are defined as: 

5 2 C(t,t') 



^ Jti,t 2 <t 
R 2 (t,t') = \ [ dhdt2 

* Jt \ .to.<t 



cos(uti) cos(ujt 2 ). (10) 



E=0 



5 2 R(t J t') 

COS[0Uti) COS{U>t2)- (11) 



5E(t x )5E(t 2 ) 



E=0 



The second order correlation function, C 2 (t, if), was introduced in the context of sp in-glasses 



3 {t| i 



31 



Neither 



TTI nor FDT holds for C 2 (t,t') and R 2 (t,if), which are explicit functions of both t and 
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t! . However, the response kernels appearing inside the above integrals are TTI. Using this 
property and developing the product of cosines in complex exponentials one finds easily that: 



C 2 (t,t') = cf\t - t>) + (e^ + *') c^\t - t') + c.c.) 

R 2 (t,t') = r^\t-t')+(e i ^ t+t '\f 1 \t-t') + c.c\. (12) 



B. Relation between non-linear susceptibility and second-order response 

In the following we aim at establishing a relation between the second order response func- 
tion defined above and the non-linear susceptibility. By definition, the electric polarization 
is given by the convolution of the response function with the external field: 

P{t) = I dt' R(t,t')E{t'). (13) 

J — oo 

Therefore, using Eqs. (jS]) and ([H]), we simply get that: 

EP 1 (t)= I dt' R {t-t')E{t'), (14) 

J — oo 



and 

-t 



EP 3 {t) = ( dt' R 2 {t,t')E{t'). (15) 

J — oo 

Thus, the component of order E 3 of the polarization (related to ^3) turns out to be just the 
convolution of the field with the function R 2 (t, t') defined in Eq. (TTTJT) . As a consequence, us- 
ing the above expression, together with Eq. ffl2]) . we find, for an oscillating field at frequency 
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P 3 {t) = 5?{e^ (f { o\u) + f{ w) (0)) + e^ff ^)} , (16) 

where we denoted f^\u') and r^\u') the semi- Fourier transform (with respect to r) at 
frequency uo' of the coefficients appearing in the Fourier expansion, 7q (t) and (r) . The 
previous equation allows us to establish a general relation between Xsi 1 ^) an d the Fourier 
transform of r^-^r): 

/"OO 

Xa (w) = 4rf } (2^) = 4 / dr e" 2 ^ rf ] \r) . (17) 

Jo 

This relation will be very useful. Using scaling arguments we will now obtain the critical 
behavior of r± (r) within MCT. The relation above will then allow us to obtain straightfor- 
wardly the scaling behaviour of Xzd^)- 
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III. LOW FREQUENCY LIMIT AND RELATIONSHIP TO DYNAMICAL LIN- 
EAR RESPONSES 



A. Low frequency limit 

In the following we focus on the evolution of the correlation and response function when 
the period of oscillation of the external field is much smaller than the relaxation time, 
in other words in the low- frequency limit ur a <C 1. By definition all degrees of freedom 
relevant for this time-sector of the response or correlation relax on timescales much smaller 
than oj~ l . As a consequence the correlation/response functions are expected to be given 
by their equilibrium expression in the presence of a quasi-constant external field Ecos(uit). 
Therefore, in this regime: 

C(M') = C eq (t - t' ', E cos{ut)) (18) 
R(t,t') = R eq (t-t',Ecos(ut)). 

Since we are interested in the small E behavior, we can expand the above expression up to 
second order in E. For the response function, for instance, this yields: 



Req(T, E COS (ujt)) « Rq(t) + 



E 2 cos 2 (ut) d 2 R eq (r,E) 



(19) 

E=0 



2 dE 2 

where Ro(r) is the unperturbed equilibrium response function, and the derivative is com- 
puted with respect to a constant external field. Comparing the last equation with Eq. f|T2|) 
in the stationary regime (t, t' — > oo with r = t — i! finite) we find a very simple expression 
for the n = 1 component of the expansion r[ (r) in the regime UT,UT a 1: 



(u)/ s 1 d 2 R eq (r,E) 
i ( r ) = o 



dE 2 

An analogous relation holds for the correlation function: 



(20) 

E=0 



1 d 2 C eq (r,E) 



dE" 1 



(21) 

E=0 



These general results provide important insights to understand the behaviour of the non- 
linear susceptibility. First of all, since the correlation and response functions appearing in 
Eq. ([TBI are defined in equilibrium in presence of a constant field, they must obey FDT. 
Therefore one can establish a sort of generalized Fluctuation-Dissipation relation between 
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the second order correlation and response functions, which reads: 

rf » (T) = _18fgM (22) 

which is however only valid in the low frequency domain ur,ujT a <C 1. 
B. Relationship with dynamical linear responses 

The results of the previous Section allow us to establish an interesting relation between 
the second order correlation and response functions and the dynamical res pon se Y t(t) = 



TdC eq (r) I dT that was recently introduced and extensively studied in [19|, |20|, |26|, |27|, |28|, |30| , 
in particular in relation with the behaviour of the four-point dynamical correlation function. 
The key idea is that in the glassy dynamics regime, the equilibrium correlation function 
C eq {r) satisfies to a good approximation the time-temperature superposition (TTS) princi- 
ple. This means that the correlation function for different temperatures, densities, external 
fields, etc., can be written as a function of r/r a (T, p, E), and the whole T, p, E dependence 
is captured by the structural relaxation time T a (T, p, E). This becomes actually an exact 
statement within the a-regime of MCT, when the system approaches the dynamical critical 
point. In this case, the dynamical critical temperature is expected to show a quadratic 
dependence on the external field (for small fields) of the form: 

T MC t{E) « T M ct{E = 0) + kE 2 . (23) 

Close to the critical point, a small field changes slightly the critical temperature. Since the 
only thing that matters for the critical behaviour is the distance from the critical point, one 
finds that applying a small field is equivalent to a small change in temperature. Note that 
this implies that the relations found below carry over, within MCT, to the /3-regime as well. 

More generally, since r a (T, p, E) is expected to be an even function of E because of 
the up-down symmetry, it should rather be written as r Q (T, p, B), with = E 2 . Then 
approximate TTS immediately leads to: 

dC eq (r) ^ drJdQ dC eq (r) 

dQ ~ dr a /dT dT ' 1 ' 

Around E = one has d 2 C eq /dE 2 = 2dC eq /dO; using the above results one finds the very 
interesting relation (valid for lot < 1): 

c { ?\r)~\ d -^l = ^XT{T\ (25) 
12 



where k = fj@-/|j!f-- Analogously one finds for the response function (using the FDT relation 

my 

rf\r) = - JL . (26, 

Using Eq. ffTTl) one directly finds the relation between % 3 and £t given in Eq. ([I]). 

One can in fact obtain a slightly more general relation provided the Fourier transform of 
Eq. (fI7|) is dominated by the region ujt <C 1. One then finds with some degree of generality 
that: 

~ dQ , (27) 
and finally, using TTS, Eq. ([2]). This result is important because it establish a firm link 
with the linear dynamical responses that are often used to evaluate dynamical correlations. 

More generally, the amplitude of the correlation function also depends on temperature 
and electric field: 

C eq (r)=A(T,E,..)c eq (-), (28) 

and the derivative of A brings extra contributions that affect the above equalities. In glassy 
systems, the relaxation time r a is usually most sensitive to external parameters, and it is 
reasonable to discard these corrections, except at zero frequencies where the above contri- 
bution is in fact zero. The case in spin-glasses is very different, because the correlation 
amplitude itself depends critically on temperature. 

IV. CRITICAL BEHAVIOUR OF THE NON-LINEAR SUSCEPTIBILITY 
WITHIN MCT: SCALING ARGUMENTS 

In this section we analyze the behaviour of the non-linear susceptibility using general 
physical and scaling arguments. These results can be confirmed by an exact analysis of the 
schematic Mode Coupling (p-spin) equations. For sake of clarity of this paper, the technical 
aspects related to the derivation of the schematic MCT equations in presence of an external 



oscillating field and their analysis will presented in a separate publication 33jj. 



A. Technical preliminaries 



The results above establish a clear connection between the dynamical response xt and 
the non-linear susceptibility %3- in the following, we will exploit the consequences of this 

13 



connection within the MCT framework, using scaling arguments. Note that we will im- 
plicitely assume, as previously done in the literature, that the results obtained within MCT 
for density correlation functions carry out to polarization fluctuations. In fact, experimen- 
tally, it has been established that the dielectric susceptibility probes the glassy dynamics 
as well as the density correlation functions, see e.g. 34j- A quantitative theory of dielectric 
polarization fluctuations and their coupling to density fluctuations in the slow dynamics 



regime would be certainly very involved due to the presence of Onsager cavity fields 



35]. 



Let us now recall the behaviour of Xt(t) within MCT |20l . |30[, where two critical re- 
laxation regimes occur close to the MCT transition: the /3-regime, with relaxation time 
Tp ~ Toe" 1 / 2 ", and the a-regime, with relaxation time r a ~ r^e" 1 / 26 ^> Tp. 

The critical properties of Xt(j) have been derived in terms of scaling functions both in 
the a and /5-regimes {^(J. Using these results, the relations established above and the FDT 
of Eq. (122]) . one can obtain the scaling behaviour of c^\t) and t^\t) close to the MCT 
transition in the different time regimes. The strategy is the following: we start with very low 
frequencies, where we know the scaling behaviour from the relation with xt, and extend it 
to the whole regimes assuming that critical scaling holds. Finally, from the scaling behavior 
of Ci(t) and r[ (t) we will obtain the one of Xst^)- We will first analyse the /3-regime and 
afterwards the a-regime. 



1. The (3 -regime 

Let us first analyze the regime lot <C 1. In this case, one can use the results (I25II26P valid 
in the regimero <C r ~ Tp <C r Q and the known behavior of Xt Q, 3(| to find: 



c 



1 



r)~ Ve C^-j (29) 

where the scaling function (x) behaves asymptotically as x a for 1 < 1 and as x b for i> 1, 
whereas rp{x) behaves for i<1 and for x 3> 1. 

On the other hand, for very large frequencies, such that lot ^> 1, the field oscillates so 
fast that the system has no time to respond and one expects vanishing (for lot — > 00) small 
second order correlation and response functions. 
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Finally, for ujt of the order of one, the most general scaling behavior in the /3-regime 
generalizing the one above reads: 

c< "' w M W^ r ) (30) 
^w-^^fe^)- (31) 

The most general assumption compatible with previous results is a factorized form for fp 
and gp in both regimes t> Tp and r <C Tp. In the late /3-regime, one has (with L for 'late'): 

cfV) ^ 4= f-V^(^) (32) 



(w)/ \ 1 1 / r x 



g L p(ur). (33) 



where both functions fp,gp tend to a constant when their argument ujt is small, and tend 
to zero when ujt is large. In the early /3-regime, a similar result holds, with a priori different 
scaling functions fp,gp (E for 'early'): 



1 / t 



a 



1 1 / \ o— 1 



/ e T /3 \Tf) 

cE „E 



-T- - (35) 



These new functions fp,gp again tend to a constant when their argument ujt is small, and 
tend to zero when ujt is large. 

In this regime, the explicit dependence with r and uj occurs only through the rescaled 
time and frequency r/rp and LOTp. In the rest of the text, we will frequently use the variables 
t — t I Tp , uj = ujTp and x = ujt. 



2. The a-regime 



The behaviour of cf\r^ in the a-regime follows similar scaling laws. When ujt <C 1, 



the results of 



20 



30| allow one to obtain: 



1 T 
-Pa — 

6 \T n 



(36) 



rj (t) ~ r c 
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The matching between the late /3-regime and the a-regime determine the asymptotic be- 
haviour of the scaling functions defined above. One finds that c a (x <C 1) behaves as x b and 
r a (x <C 1) as x , whereas both functions tend exponentially fast to zero for i>1. 
When lot is not small, the scaling behaviour in the a-regime reads: 

4">M*iL(J-) (37) 

rS"»W « --gj-,ur). (38) 
In the early a-regime, that is when r <C T a , one finds: 

6 ^f(^). (40) 

Furthermore, by requiring the matching between the two regimes of large t/t^ and small 
t/tq, one finds that the scaling functions /<f (x) and g^ix) are the same as /J'(x) and gfj(x). 

As in the /^-regime, the explicit dependence with r and u; occurs only through the rescaled 
time and frequency r/r a and ujr a . In the rest of the text, we will frequently use the variables 
t = r/r a , To = uJT a (and also x = lot). 

In Fig. [2] we show a sketch of the behavior of (c^ (t)| that summarizes all our previous 
findings. 

B. Scaling behavior of X3( w ) 

In the previous section we have determined the scaling forms governing the time and 
temperature dependence of r[ (r) and c{\t). Using these results we can now easily analyze 
the critical behaviour of the non-linear susceptibility, by computing its Fourier transform at 
frequency 2lo, according to Eq. ( 1TTI) . We again focus in turn on the /3-regime and then on 
the a-regime, before commenting on the zero and infinite frequency limits. 

1. The (3 -regime 

Let us first consider probing frequencies of the order of the inverse of the /3-relaxation 
time. We set Co = lot^ and assume that only the /3-regime of ?i (t) contributes to Eq. (JTTj). 

1(3 



1 


1 




€ 








FIG. 2: Sketch of log (c^ (r)| as a function of logr. In the u> — ► limit Cj (r), behaves as xt(t"), 
i.e., it scales as 1/v^ { T / T p) a hi the early /3-regime and as l/y^ (t/t^) 6 in the late /3-regime (or, 
equivalently, as l/e(r/r Q ,) f ' in the early a-regime). In the a-regime \c[ (t)| reaches a maximum of 
order 1/e. At finite frequency oj, \c\ (t)| drops to values of 0(1) as r > 1/cj (red curves). 

This assumption can be fully justified by the exact analysis of the schematic Mode Coupling 



equations 



Indeed one can show that r^\r ~ r a ) is small due to the fact that the 
scaling function g^(x), introduced in Eq. (1401) , vanishes as l/x l+b at large x. Therefore, for 
probing frequencies of the order of the inverse of the /3-relaxation time, the contribution to 
the non-linear susceptibility coming from the time integral in the a-regime is negligible in 

Eq. jrr 



One then finds that: 

xM^L^pH*) (4i) 

v e 
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where the function J-{x) is defined as: 




The asymptotic behaviour of the scaling function JF can easily be obtained from the results 
of the previous section. One finds: 

/>oo 

T(£o) ~ ACo- b duu^e-^g^u) ujt p < 1 (42) 
Jo 

POO 

~ 4iu~ a d««°- 1 e- 2<B y|(M) wr^»l. (43) 

Using the asymptotic properties of g^(u), g§(u) and the fact that a,b are between zero 
and one insures the convergence of the integrals appearing in the above equation, at both 
small and large u. This confirms that the scaling behaviour of Xsi^) m this region is 
indeed dominated by the /3-regime of r[ (t). Note that in the high frequency region the e 
dependence of Xzi^) drops out, as it should in order to match the non critical Tq 1 frequency 
regime. 



2. The a-regime 

We now consider the a-regime, where we set uj = ujr a , and again assume that only the 
scaling form of r± (r) in this same regime, Eq. (j38j) . contributes significantly to Eq. (fTTj) . 
Indeed, the contribution due to the time integral in the /3-regime is at least a factor y/e 
smaller than the one coming from the a-regime, and yields a subleading contribution to the 
critical behavior of Xsi^)- We then find that the non-linear susceptibility scales as: 

xM-\g&) (44) 

where the function Q{x) is defined as: 

Q(P) = = &ue-^-g E a (^u) 
to J \u J 

Clearly, because of the matching of ~g a for small first arguments with gp at large first argu- 
ments, we find that the scaling of the early a-regime (ujr a ^> 1) of the non-linear suscepti- 
bility matches with that of the late /5-regime (corp Cl), with: 

X s(co) oc e ^' 2a u- h . (45) 
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3. Low frequency limit 



In the low frequency limit, one finds that X3( w ) decreases from its peak value e _1 reached 
for uo ~ l/2r Q to a non critical, finite value given by Eq. (J27J). As discussed in contrary 
to the case of spin-glasses, the non-linear susceptibility is critical only for small but non zero 
values of the frequencies. Zero frequency corresponds to a static equilibrium response (or 
correlation, via FDT). In glasses, these are not expected to have any critical behavior. 

In particular, in the low frequency limit, one can expand Eq. (fTTj) up to second order in 
lot. Using Eq. (TT9]) we have: 

k r dr (1 - 2iur - 2u 2 r 2 ) (46) 
Jo oT 

« «^{Xi(0) (l-2iLuA 1 T a -4uj 2 A 2 T 2 )}, 

where the zero frequency limit of the linear susceptibility, Xi(0), equals the static polariza- 
tion fluctuations (along the z-axis) divided by temperature, N(P 2 )/T. A\ and A 2 are two 
temperature-independent constants defined as A\ = / °° dss* r eg (s) (here we have used again 
the time-temperature superposition principle, writing R eq {r) = Xi(0)r - eq (T / 'r a ) / V a ) . The last 
equations allows us to determine the low frequency behavior of the real and imaginary part 
of the non-linear susceptibility: 

9f(X3H) « J B 2 cu + 0((^r a ) 3 ), 
with B 1 = -4KA 2 Lj 2 d(xi(0)T 2 )/dT > and B 2 = -2kA x uj d(xi(0)r a )/dT > 0. 



^. Large frequency limit 



At very large frequencies (very small timescales) the non-linear susceptibility is vanishing 
because the system has not enough time to respond to the oscillating field. One could argue 
that the analysis of Eq. ([6]) at very large frequencies yields: 



Ps(t) 



5 3 P 



5E 3 (0) 



E=0 



dti (e 



iasti _|_ iu>t\ N 



L-'O 



As a result, at very large frequency the non-linear susceptibility behaves as: 

1 S 3 P 



oo 



(47) 



(48) 



E=0 
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This analysis is oversimplified and assumes analytic properties of the function sE^SEfo^Efo) 
that are not granted and may depend strongly on the microscopic dynamics. For instance, 
in the case of Ising spins with a Monte Carlo heath bath dynamics one can easily verify 
that the previous arguments do not apply and the large frequency behavior is proportional 
to The conclusion is that the high frequency behavior depends on the underlying 

microscopic dynamics and, likely, on the type of the non-linear response considered. In the 
case of non-linear dielectric susceptibility the underlying microscopic dynamics should be 
provided by Langevin equations for dipoles in a non-polar solvent (this is an approximation 
since at extremely high frequency inertia effects will play a role). To work out the high 
frequency behavior one can neglect interactions with other dipoles and the coupling to 
structural relaxation. Thus, the analysis of the non-linear response of a single dipole in a 



non-polar solvent worked out in 
discussed above. 



361 ] should apply. The outcome is the l/(iu) 3 behavior 



C. Beyond schematic Mode Coupling: general considerations and non-linear re- 
sponses to time and space inhomogeneous fields 

The previous scaling arguments are rather general. We emphasis here that they can 
be fully derived from a rigourous analysis of schematic Mode Coupling equations, in the 
framework of the spherical p-spin model in presence of an oscillating external field. Analysing 
the equations up to second order in the external field, it is possible to determine the critical 
behaviour obtained using scaling arguments in the previous sections. 

Beyond the schematic approach, one may wonder about specific but important details 
such as the role of conserved variables like energy or density on the above results. As 
recalled in the introduction, we know that these conserved variables can dramatically change 
the scaling behaviour of X4 f° r example, which also diverges as (T — T c ) _1 within a p-spin 
framework with Langevin dynamics, but diverges as (T — T c )~ 2 when the contribution of 



conserved variables is taken into account 
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20| . From a diagrammatic point of view, this 



is due to the presence of 'squared ladder' diagrams which gives the dominant contribution 
to X4- One can check that due to the causality of the response functions, these diagrams in 
fact are absent when one computes the non-linear susceptibility and the above results are 
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expected to hold for a bona fide MCT theory of liquids. Thus, the result [l|: 

Xsi^-e^GiuTa), (49) 

is expected to hold even beyond schematic MCT (£ is the dynamical correlation length which 
also appears in xt)- A complete proof could be obtained generalizing the Inhomogeneous 
MCT calculation of |30J to account for a space and time dependent source term, that would 
describe the non-linear response to an oscillating field with wave-vector q and frequency u. 



In a first attempt, we adapt the previous scaling arguments to the explicit results of [30j on 
the wavevector-dependent dynamical response. This allow us to obtain the critical behavior 
of Xzi^-, <?)• I n the /3-regime one finds: 

Xs(w, q) = £ 2 Hp{ujrp, q£) £ = e~ 1/4 , Tp = 

where the scaling function Hp(x, y) is equal to T{x) for y — 0, i.e. for a uniform electric field 
one finds back Eq. (I4ip . For large y one expects a power law behavior such as h ^f l (where 



hp{x) is a certain scaling function). As discussed in 30] this is needed to cancel out the 
diverging prefactor £ 2 and match the critical behavior to the non-critical one taking place 
for q oc 0(1). The asymptotic behavior with respect to x is identical to the one already 
described for homogeneous fields. Very small x correspond to the matching between a and 
j3 regimes. Since in the a regime X3 is expected to diverge as £ 4 , the matching imposes the 
behavior at small x: (where hp(y) is another scaling function). For large x values, the 

field varies so rapidly that the system has not enough time to adjust and to respond to the 
field. Again, in order to cancel the diverging prefactor and match the non-critical behavior 
one expects a large x behavior such as hp Jf (where hM{y) is a third scaling function). It 
would be interesting to specify in more details the shape of the scaling functions hp, hp, hp. 
In the a regime one expects: 

q) = £ A H a {uT a , q£) £ = e~ 1/4 , r a = e™^"^ 

where the scaling function H a (x, y) is equal to G(x) for y — 0, i.e. for a uniform electric field 
one finds back eq. (jSJ). The same kind of arguments used above suggests for large y a power 
law behavior: H a (x,y) ~ ha ^ (where h a (x) is a scaling function). For very small x the 
scaling function vanishes in order to match the x = value corresponds to the non-critical 
(non diverging) static non-linear susceptibility. For large x values in order to match the j3 
regime one expects a behaviour such as ha ^ (where h^(x) is another scaling function). 
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At this point, one may debate about the validity of the power law divergence of the relax- 
ation time described by Mode Coupling Theory and of TTS. However, it is often observed 
that there are actually regimes, when the dynamics starts slowing down strongly, which are 
well described by MCT (and thus TTS). In these regimes we expect our predictions for ^3 
to hold, as it seems to be the case for % 4 [3?}]. Furthermore, it was recently shown that 
MCT may be corrected including higher order terms. It was found that these corrections 
only affect the values of exponents a and b, and T c , but do not affect strongly scaling func- 



tions 
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391 ] . This suggest that the MCT mechanism for describing the slowing down is 



rather robust and that the MCT regime could in principle being expanded at the price of 



changing the exponents (in particular the one coverning the relaxation time, see e.g. 39() 



V. CONCLUSION 

In this work, we have studied in detail the non-linear response of supercooled liquids. 
Although we are able to provide precise statements within a Mode-Coupling approach, some 
of our results are in fact more general and only require Time- Temperature Superposition to 
hold. In particular we expect our results to hold in generalization of Mode Coupling Theory 



such as 



40]. 



An important theoretical result is the relation (T2T1) between the non-linear response X3(cu) 
and the temperature derivative of the usual linear susceptibility, dxi{2u) / dT , valid at small 
frequencies (smaller than the inverse relaxation time). This bridges the gap between non- 
linear responses and probes of dynamic correlations such as three- and four-point correlations 
and temperature (or density) derivative of standard two-body correlations and response: 
they are all different facets of the same underlying physical phenomenon. 

For larger frequencies - of the order of the inverse of the relaxation time - we have obtained 
a complete set of scaling predictions concerning the critical behavior of Xsi 1 ^) within MCT. 
The main results are summarized in Fig. [TJ Five different frequency regimes are identified: 
Xs(u;) exhibits a peak around frequencies of the order of half the inverse of the structural 
relaxation time of the system r a . The height of the peak grows as (T — TJ^ 1 (or equivalently 
as £ 4 ) as the critical temperature is approached. For higher frequencies, Xsi 00 ) decays as 
power laws, with an exponent equal to —b in the late /3-regime, to —a in the early /3-regime, 
and finally to —3 at high frequencies. 
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Our results should be directly applicable to the non-linear dielectric constant of molecular 
glasses in the weakly supercooled regime where MCT is expected to be relevant, and for 
describing the non-linear compressibility or more general non-linear rheological responses of 
hard-sphere colloids close to the glass transition, where MCT does a fair job at describing 
their relaxation properties. However, it is well known that MCT fails for deeply supercooled 
liquids, when activated events start playing a major role in the relaxation. The detailed 
shape of x^{ u ) would require a full theoretical description of the dynamics in this regime, 
which is unavailable to date. Still, the general low frequency relation between % 3 and dxi/dT, 
supplemented with the property of Time- Temperature superposition, allows one to give a 
firmer basis to the scaling relation conjectured in namely that: 

XsM ~ xlG(^T a ), (50) 

where Q is a scaling function, and x% ^ d\n. T n f dh\T is the peak value of the temperature 



derivative of as measured in [26|, |27|, [28(. Following [19|, |20|, [26], we expect Xt to 



increase as a power of the dynamical correlation length £ (T) . The detailed shape of Q would 
obviously be worth knowing in order to compare with upcoming experimental results. As a 
guide, we give the result obtained assuming a Havriliak-Nagami form for the susceptibility 
and the validity of the relation between X3 an d dxi/dT at all frequencies, which has no 
justification apart from suggesting possible fitting functions. One finds: 

Shn(u) = {1 + {m)by+c , (51) 

where b, c are fitting exponents. 

Among open problems worth investing is the extension of the present theory to the aging 
regime of glasses and spin-glasses. From an experimental point of view, a detailed study 
of the role of the electric field on the glass properties of dipolar liquids (such as glycerol) 
would be very interesting. For example, the evolution of the glass transition temperature 
as a function of the field E would allow one to measure the proportionality coefficient k 



appearing in Eq. (1271) . In spin-glasses, a detailed measurement of Xz{^) would allow to 
shed light on the existence of spin-glass transition at non zero field, as argued in [l]. The 
behaviour of x 3 (c t ;,t u ,) in the aging phase would furthermore be a very useful probe of the 
aging process in spin-glasses, in particular during rejuvenation cycles. Numerical simulations 



of X3(ui,t w ), using the zero-field techniques developed in [3l|, |32j, would be worth pursuing. 
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Figure 1 Sketch of log|x3(w)| as a function of logo;, showing five different frequency 
regimes: ur a <C 1, ur a ~ 1, Tp/r a <C urp <C 1, UTp ^> 1, ujtq ~ 1. Note that the low 
frequency limit is non zero but much smaller than the peak value for T close to T c . 

Figure 2 Sketch of log|c^(r)| as a function of logr. In the uj — > limit (r), 
behaves as Xt(t), i.e., it scales as 1/y/e (t/t^Y in the early /5-regime and as 1/y/e (t/t^) 6 
in the late /3-regime (or, equivalently, as l/e(r/r a ) 6 in the early a-regime). In the a-regime 
|c^(r)| reaches a maximum of order 1/e. At finite frequency u, \c\ (t)| drops to values of 
0(1) as r > \/uj (red curves). 
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